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Abstract 

Within the program of holographic renormalization, we discuss the 
computation of three-point correlation functions along RG flows. We 
illustrate the procedure in two simple cases. In an RG flow to the 
Coulomb branch of J\f = 4 SYM theory we derive a compact and 
finite expression for the three-point function of lowest CPO's dual to 
inert scalars. In the GPPZ flow, that captures some features of M = 1 
SYM theory, we compute the three-point function with insertion of two 
inert scalars and one active scalar that mixes with the stress tensor. By 
amputating the external legs at the mass poles we extract the trilinear 
coupling of the corresponding superglueballs. Finally we outline the 
procedure for computing three-point functions with insertions of the 
stress tensor as well as of (broken) R-symmetry currents. 



1 Introduction 

Holographic renormalization [1, 2] is a powerful tool that produces finite 
correlation functions and allows explicit checks of (anomalous) Ward iden- 
tities along RG flows dual to asymptotically AdS domain walls. One of the 
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main achievement of the program is an efficient and consistent algorithm for 
computing finite one-point functions in the presence of external sources. By 
sources we mean the coefficient functions, two per each bulk field compo- 
nent, which are not fixed by the near-boundary analysis [3, 4]. One of those 
corresponds to an operator deformation of the fixed point theory, while the 
other corresponds to giving a VEV to an operator [5]. 

The non-local and non-linear relation between VEV and source which 
amounts to requiring regularity of the classical solution of the non-polynomial 
bulk field equations is hard to find in general. Luckily in order to compute 
two-point correlation functions it is enough to know this relation at the lin- 
ear level, e.g. by solving the linearized fluctuation equations around a given 
domain wall solution. Disentangling the mixing between longitudinal/trace 
modes of the stress tensor and the active scalar (s) proved subtler [1, 6], 
than initially expected [7, 8]. Similar difficulties appearing in the mixing of 
the broken symmetry currents with Goldstone bosons were disposed of in 
[2, 9, 10]. Transverse modes are easier to analyze since their fluctuations can 
be expressed in terms of an auxiliary massless scalar [11, 12]. In particu- 
lar, after getting rid of non-canonical kinetic terms, (transverse) vector fields 
display a universal effective mass M 2 = —2A" [11, 13, 14] that is manifestly 
positive thanks to the convessity of the scale factor A(r) which is at the heart 
of the holographic c-theorem [15, 16]. Interesting mass spectra emerge from 
momentum analysis of holographic two-point functions but so far tradition 
has prevented the study of even the simplest interactions. 

There is little hope of computing higher-point correlation functions except 
possibly for three-point functions 3 . Higher point functions require the knowl- 
edge of the bulk-to-bulk propagator that appears in exchange diagrams. For 
three-point functions, the relevant Witten diagrams at tree-level only involve 
bulk-to-boundary propagators that are known in many cases, e.g. solving the 
linearized fluctuation equations with prescribed boundary conditions. Bulk- 
to-bulk propagators would only appear in loop diagrams that are suppressed 
by the 5-D Newton constant i.e. by inverse powers of N in the spirit of the 
AdS/CFT correspondence [17, 18, 19] 4 . Thus restricting to the planar limit 
captured, at strong 't Hooft coupling, by the supergravity approximation one 

3 Following the procedure of holographic rcnormalization, three-point functions at the 
conformal point dual to pure AdS have been computed by D. Freedman, U. Gursoy and 
K. Skenderis as reported in [23]. 

4 For recent reviews see e.g. [20, 21, 22, 23] and references therein. 
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can go a step beyond the tradition and address the problem of computing 
three-point functions. By 'computing' we actually mean writing down closed- 
form expressions that allow one to extract interesting information about the 
interactions in particular regimes. The only ingredient that is missing from 
previous analyses is the exact form of some bulk cubic couplings along the 
flows and of some additional counterterms. In this note we fill this gap and 
compute simple three-point functions of both active and inert scalars in the 
two most studied holographic RG flows: the GPPZ flow [24] and the CB flow 
[29, 12]. 

In order to fix the notation we briefly review the procedure of holographic 
renormalization in Sect. 2. We then describe the general strategy for com- 
puting three-point functions in Sect. 3 and collect some basic formulae for 
the two flows in Sect. 4. The only real novelty that deserves a special men- 
tion is the identification of additional logarithmically divergent counterterms 
in the GPPZ flow when the bulk field dual to the gluino bilinear comes into 
play. 

In Sect. 5 we specialize our analysis to the CB flow. With little effort we 
write down a compact expression for the three-point function of inert scalars 
belonging to the (3,3) representation of the 5*0(4) subgroup of SO (6) pre- 
served along the flow. We check that our result is finite for generic momenta 
and try to convince the reader that it displays the correct UV behavior for 
large momenta. For comparison and in order to give a flavor of the com- 
plexity of the problem, we compute the Fourier transform of the conformal 
three-point functions of scalar primary operators which is not easy to find in 
the existing literature [25]. We thus report the relevant formula for A = 2 
in an appendix. 

In Sect. 6 we analyze the three-point function with one insertion of 
the active scalar that interpolates for the trace of the stress tensor and two 
insertions of the inert scalars, which is a "caricature" of the gluino bilinear 
in the GPPZ flow. As expected the resolution of the mixing with the stress 
tensor at the required order in the fluctuations turns out to be quite laborious 
but straightforward. The naive result is logarithmically divergent in the 
contact terms that are cancelled by the additional counterterms previously 
identified. By amputating the external legs we extract the finite and rather 
explicit trilinear coupling of the relevant superglueball states. 

Ideally one would like to compute three-point correlation functions with 
insertions of the stress tensor and/or of (broken) R-symmetry currents. We 
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plan to return to this problem soon. In the concluding Section we outline 
the procedure for computing such three-point functions and speculate how 
similar computations, following or extending when needed the approach of 
holographic renormalization along the lines of [14], should be feasible for the- 
ories which are dual to full-fledged string solutions possibly with logarithmic 
corrections even in the UV [26, 27, 28]. 



2 Holographic Renormalization strikes back 

In this section we fix our notation. Since there only minor changes with 
respect to [1, 2, 23], readers familiar with holographic renormalization may 
skip directly to Section 3. 

To be specific we focus on asymptotically AdS domain walls which are 
dual to RG flows in d — 4. The set of fields driving the flow together with 
some 'spectator' fields represents a consistent truncation of some gauged 
supergravity in D = 5, which in turn is the low-energy effective theory of 
some superstring compactifications on an Einstein manifold with stabilizing 
fluxes. As usual, we consider Poincare invariant solutions of the form [15, 
29, 16, 24, 12] 

ds 2 = e 2A(r) dx-dx + dr 2 $ = $(r) (2.1) 

where <3> denote one or more canonical scalar fields that drive the flow. The 
AdS boundary (field theory UV) is at r — > +oo. The (IR) interior presents 
a naked singularity unless the flow terminates at another (super) conformal 
fixed point. Preserving some supersymmetry reduces the problem to solving 
first order 'gradient flow' equations 

A! {f) = -f &(r) = 9 -d^W (2.2) 

3 2 

where g = 2/L is the bulk gauge coupling and L is the characteristic AdS 
scale. The scalar potential reads 



q 2 

V = — 



2 I 0$ j 3 



(2.3) 



To extract some physics from the domain wall solution one may not forgo 
a near-boundary analysis of the non-linear field equations [4] . To this end it is 
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customary to introduce a different radial variable z such that dr = —Ldz/z, 
i.e. z = exp(— r/L). The scale factor behaves as A(r) m r / L m — log z near 
the UV boundary at z — 0. When exp(2A) exp(+2r/L) as r — > — oo one 
has an AdS 'horizon' with L < L that describes a different fixed point at 
the endpoint of the flow [29, 30]. Barring string corrections or uplift to D — 
10 [31, 32], many interesting solutions display a naked singularity at some 
finite value of the radial variable. Depending on the behavior of the active 
scalar field(s) near the boundary one has flows along which (super) conformal 
invariance is either spontaneously broken, by the presence of operator VEV's, 
or explicitly broken, by a relevant deformation of the fixed point action. 

For scalar fields dual to operators with UV dimension A, i.e. with AdS 
mass (ML) 2 = A(A — 4), the near boundary expansion, that reads 

$(z,x) = z 4 - A ^ (0) (x) + z(j> {1) {x) + ...) + z A (4> (0) (x) + z4>(i)(x) + ...) (2.4) 

displays two a priori independent asymptotic behaviors. The iterative solu- 
tion of the field equations allows one to express the higher order coefficient 
functions 0( n ) and 0( n ) in terms of 0(o) and (f)my For A an integer, only co- 
efficient functions with 'even' indices appear and it is convenient to change 
radial variable to p = z 2 . Moreover the two dominant behaviors differ by 
an integer and additional logarithmic terms are needed in order to satisfy 
the field equations 5 . Logarithmic behaviors can be associated to conformal 
anomalies that, being local, can be written in terms of (f>(p), playing the role 
of an operator source, and a finite number of derivatives thereof. It is not 
difficult to see that 0(o) is related to the VEV of the corresponding operator 
[5, 4]. 

Naive extension of the AdS/CFT correspondence to holographic RG flows 
suggests one to compute the on-shell action as a functional of the boundary 
data and then take functional derivatives in order to compute correlation 
functions. However the unregulated on-shell action is infinite due to the 
infinite volume of the bulk. A convenient regularization consists in a radial 
cutoff p > e that allows one to easily isolate a finite number of terms diverging 
as e — > 0. For the case of scalars coupled to gravity the on-shell action takes 

5 In the very special case A = 2, even the dominant term contains a logarithm <&(/?, x) = 
p\ogp<p (Q) (x) +p(f>( Q) (x) + .... 
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the schematic form 



>W0(o),£(o);e] = / d 4 x v ^y + ^ + ... - loge a {2u) + O(e ) 
J p=e Lee 

(2.5) 

where v = 2 and a^k) are local functions of the sources, 0(o) and g^o), that 
appear in the expansions of $ and of the metric 

9(o)ij = 9(o)ij + P9(2)ij + P 2 9(4)ij + P 2 logp/i(4)ij + ••• (2.6) 

The counterterm action S^t^, 7; e], chosen in such a way as to cancel the 
divergent terms in S Teg [(j)(o), g^y, e], is to be expressed in terms of the fields 
<3>(x, e) 'living' at the regulating surface p = e and of the induced metric, 
7ij = gij(x, e)/e. This is required for covariance and follows from a straight- 
forward but often laborious "inversion" of the expansions (2. 6), (2.1). The 
counterterms are universal in that they make the on-shell action finite for 
any solution of the bulk field equations with given boundary data. The coun- 
terterms are however different for different (consistent) truncations, e.g. for 
different potentials V($). So far we have described a "minimal" scheme in 
which only the divergences of S Teg are subtracted. As in standard quan- 
tum field theory, one has the freedom to add finite invariant counterterms. 
These correspond to a change of scheme aimed for instance to restoring some 
(super) symmetry. 

For later purposes, it proves convenient to define a subtracted action at 
the cutoff 

Ssub[$,7;e] =S reg [0(o),0(o);e] + S , c t[$,7;e]. (2.7) 

The subtracted action has a finite limit as e — > 0, and the renormalized action 
is a functional of the sources defined by this limit, i.e. 

Srcn[(f>(o),9(o)] = lim5' sub [$,7;e] (2.8) 

The distinction between S su i, and S rcn is needed because the variations re- 
quired to obtain correlation functions are performed before the limit e — > 
is taken. In particular, the expectation value of a scalar operator in the 
presence of sources, defined by 

(O,) = (2.9) 
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where = det(g^o)ij) can be computed by rewriting it in terms of the fields 
living at the regulating surface 6 

«w - as (i^fe) (2 - 10) 

where 7 = det^j). In general one can prove that 

(Cg = (2A - 4)0 ( o) + local (2.11) 

where the local terms are completely fixed by the choice of 0( O ). 

The hard problem is to determine the non-local relation between the 
operator VEV 0(o) and the operator source 0(o) that guarantees regularity of 
the solution in the deep interior (IR regime) where the supergravity solution 
may develop a singularity. Traditionally only two-point function have been 
explicitly computed along holographic RG flows [1, 2, 12, 9, 6, 10, 30]. These 
only require the non-local relation at linear order 

^) = i^zri J d4x>c ^ x - (2-i2) 

or better, in momentum space, 

fe)(P) = ^4^(p)0 ( o)(p) (2.13) 

that is established by solving the linearized fluctuation equations with pre- 
scribed boundary conditions, i.e. regularity in the interior. Upon differenti- 
ation one gets the two-point correlation function 

{0^)0 ^{x')) = C 2 (x - x') + contact terms (2.14) 

3 Three-point functions: general strategy 

Our aim in this section is to sketch the further necessary steps that one 
should take in order to compute three-point functions. For simplicity we 
detail the procedure for inert scalar fields with cubic couplings, a case that 

6 For scalars of UV dimension A = 2, an additional loge is needed in this formula. 
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finds application in Section 5. Mixing with metric fluctuations or with vector 
fields is not discussed in this section. We will address this issue in Section 6. 

For three-point functions one needs to take into account quadratic terms 
in the field equations. To this end consider the Euclidean action for a real 
inert scalar around a domain wall 7 



N 2 r 1 11 

S = — 2 J d'xdpVG[-Z{p)m 2 + + 3 T (P) $3 + -] (3-15) 



The effective mass M(p) and the cubic coupling T(p) may depend on the 
radial variable p via their functional dependence on the active scalar (s). The 
kinetic term Z(p) may also depend on p. In the cases we explicitly discuss 
here scalar fields are canonical. In order to keep our formulae as compact as 
possible we thus assume Z — 1 henceforth. 
The non-linear field equation is of the form 



At the classical level the solution can be found by summing trees. Loops 
are suppressed by powers of 1 / N 2 and we neglect them since they can only be 
consistently taken into account after embedding the model in string theory. 
For a given boundary source j(x) we find 



where D(x, p;x' , p') is the bulk-to-bulk propagator, K(x,p;x') is the bulk- 
to-boundary propagator and ... denote higher-order terms in j(x). 

The linearized approximation, which is enough to compute two-point 
functions, yields 



(V^-M 2 )$ = T$ 2 + ... 



(3.16) 




(3.17) 



= K(p,p)j(p) 

in momentum space. The near boundary anlysis then gives' 



,8 



(3.18) 




(3.19) 



7 The overall constant is determined by careful reduction from D = 10 [19]. 
8 The second index denotes the order in j. 
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and 

0(o|i)(p) = 2£=lC 2 {p)3ip) (3.20) 

To quadratic order in j 

*Q)(P,P) = J dp'VGD(p,p;p')T(p') J d 4 qK(q,p')j(q)K(p-q,p')j(p-q) 

(3.21) 

Expanding both sides in powers of p yields 

0(o|2)(p) = O (3.22) 

and 

^(o|2)(p) = 2A - 4 / d P^ K (P'P) T (p) J A^(?, p)j(q)K(q - p, p)j(p - q) 

(3.23) 

since near the AdS boundary, where p ~ 0, 

£>(x, p; a/, p') « J^Kix; P) (3-24) 

Plugging this back into the on-shell action or better into the one-point 
function given above, differentiating w.r.t. j and then setting j to zero, one 
finally gets 

(0^)0^)0^)} = {2n)H(£pi) I dpjG(p)T(p)I[K(p t ,p) 

(3.25) 

modulo contact terms that are polynomial in the momenta and cancel po- 
tential divergences arising when the regulator of the radial integration is 
removed. 

In all cases we consider, the linearized fluctuation equations reduce to hy- 
pergeometric equations in a suitable radial variable w = w(p). The solution 
which is regular in the deep IR interior (which is a naked singularity of the 
'good' kind [33]) is of the form 

4> p (w, x) = w a (l — w) l3 F(a, b; c; w) exp(ipx) (3.26) 

Analytically continuing the hypergeometric function from w — 0, that rep- 
resents the deep IR interior (singularity), to w — 1, that represents the UV 
boundary, determines the bulk to boundary propagators. The reader may 
find some useful formulae in appendix B. 
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4 Our favorite holographic RG flows 



For completeness we gather in this section some relevant formulae about two 
particularly simple yet interesting asymptotically AdS domain wall solutions: 
the CB flow [29, 12] and the GPPZ flow [24]. The former is dual to an RG 
flow from M = 4 SYM at its maximally superconformal point to a locus in the 
Coulomb branch where superconformal symmetry is broken spontaneously. 
The latter is dual to an RG flow to a 'confining' Af = 1 SYM theory where 
superconformal invariance is broken explicitly by turning on relevant (mass) 
operators for the three chiral multiplets. 



4.1 Coulomb branch flow 

The first case we consider is the CB flow with n = 2, in the notation of [29]. 
The supersymmetric domain wall solution describes a continuous distribution 
of D3-branes on a disk. The relevant superpotential is given by 

W($) = -e' 2 *^ - -e 4 ^. (4.27) 

Using (2.3), one can easily compute the potential. Near $ = it can be 
expanded as 

V(*) = ~~ + ^7i $3 + °(^) ( 4 - 28 ) 

which exhibits a tachyonic mass M 2 L 2 = —4, as appropriate for a field dual 
to a chiral primary operator (CPO) of UV dimension A = 2. 

It is very convenient to express the domain-wall solution in terms of a 
new radial variable v. The solution and the relation between v and r are 

The boundary is at v — 1 and the solution has a curvature singularity at 
v — 0. The parameter £ is the radius of the disk of branes [29, 30] in the 
10-dimensional "lift" of this solution and sets the scale of a mass gap in the 
spectrum of excitations to m = £/L 2 . 
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The scalar field vanishes at the boundary at the rate 



$ » ~(1 - v ) = ~ exp(-2r/L) (4.30) 

that corresponds to the dual operator of UV dimension A = 2 taking a 
VEV. The endpoint of the flow describes M = 4 SYM theory at a locus in 
the Coulomb branch where the 5*0(6) R-symmetry is broken to 5*0(2) x 
50(4). For any finite N, the 50(2) symmetry is actually broken [2, 9]. 
The scalar that gets a VEV is the neutral singlet component of the CPO 
Qj , = Tr {X i X i )\ w that decomposes according to 20' -> (1, 1) ©(1, l)+ 2 © 
(1, 1)_ 2 © (2, 2)+i © (2, 2)_i © (3, 3) under 50(4) x 50(2). 



4.2 GPPZ flow 

The GPPZ flow [24] corresponds to adding an operator of dimension A = 3 
to the Lagrangian, namely the top component of the superpotential AW = 
E/Zf, that gives a common mass to the three M = 1 chiral multiplets 
Zj appearing in the decomposition of the M = 4 vector multiplet. The 
solution was proposed as the holographic dual of pure N = 1 SYM theory 
[24] . Although it does not capture all of the expected properties of the field 
theory, it is particularly simple and still displays some interesting features 
including a discrete spectrum of superglueballs. 

The active scalar is a singlet under an 50(3) subgroup of 50(6). A 
consistent truncation to the 50(3) singlets yields Af = 2 gauged supergravity 
coupled to two hypermultiplets describing a G 2 (2)/50(4) coset [31, 11]. After 
lengthy calculation one gets the 5d superpotential that reads 



W(<&,E) = ~ 



/2$\ 

cosh I -j= 1 + cosh(2S) 



|-i($ 2 + S 2 ) + ... (4.31) 



Near $ = E = 0, the potential admits the expansion 

V($, S) = — — — ^($ 2 + S 2 ) - ^($ 4 - 3S 4 + 6$ 2 S 2 ) + ... (4.32) 

The mass of both $ and X is M 2 L 2 = —3, indicating that the dual scalar 
operators have UV dimension A = 3. For simplicity we only consider the 
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domain- wall solution with E = 9 , and 



,2A 



U 



(4.33) 



where w = 1 — exp(— 2r/L). The boundary is at w = 1 and the solution has a 
naked singularity of 'good' type at u — 0. Since $ v / 3exp(— r/L) near the 
boundary, we are dealing with an operator deformation rather than a VEV. 
Along the flow 



W(u) = 
W H (u) -- 



3_ 

2 — u 

u 



W^u) = -Vs 



u 



u 



4 

\/3 u 



(4.34) 
(4.35) 



For our later purposes, we need the effective mass of the field £ along the 
GPPZ flow 



3 

Z 2 



2 cosh 



2$ N 



_3_ 4-3u 



u 



(4.36) 



and its cubic coupling to the active scalar $ 



d^Ml(u) = V rh 



4^3 . , /2$ N 
sinn —= 

VV3, 



8^ V / T^M 



L 2 



L 2 



(4.37) 



also along the flow. 



5 Three-point functions in CB flow 

In order to simplify our lives we only consider inert scalars of UV dimension 
A = 2, i.e. those among the 20 tachyonic scalars which are not singlet under 
the 5*0(4) x SO (2) preserved in the CB flow. In particular we concentrate 
our attention on the (3, 3)o components Q %L%R and aim to compute 

G^ R ' jUR ' kLkR (x 1 ,x 2 , x 3 ) = (Q iLiR (x 1 )Q^ R (x 2 )Q kLkR (x 3 )) (5.38) 

9 More general solutions with a non-trivial profile for the bulk field S dual to the gaugino 
bilinear have been considered by GPPZ. 
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The rationale behind this choice is that there is a unique SO (A) x 5*0(2) in- 
variant trilinear coupling among the Q's. It is proportional to £i L j L k L £i R j R k R - 
In order to avoid additional derivative couplings, it is necessary to identify 
three components of Q %LlR with canonical kinetic term around the CB do- 
main wall solution. It is well known that the icosaplet parametrizes a coset 
SL(6)/ S0(6). Taking as coset representative a symmetric matrix S with 
det(iS) = 1, the scalar potential reads 



V = - 9 -[{triS)f-2triS*)} 

S can be diagonalized by an orthogonal SO (6) transformation 

S = R{6)S{(3)R T {6) 



(5.39) 



(5.40) 



where S((3) = diag(exp(2(3 1 ), exp(2/3 6 )) with J2iPi = 
then be expressed in terms of five independent fields, 
that leads to canonically normalized scalars is 



0. The six /3's should 
A convenient choice 



Pi 
ft 

ft 
06 



x/2 



-$ + V3(X + a + uj)} 
-$ + -a-uj)\ 
-$ + y/3(-X + a - u)] 
-$ + V3(-X -a + u>)] 



2$ + VEu] 



2$ - VEv] 



(5.41) 



Taking $ as active scalar, one can convince oneself that v transforms un- 
der the lower SO{2) subgroup while A, a, uj only transform under the upper 
SO (A). The former is part of the charged singlet (1, 1)±. The latter three 
are part of the neutral ennaplet Q %h%R e (3, 3) . 
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We can proceed by switching off all the other inert scalars, setting $ to 
its background value and keeping only terms up to cubic order in A, a, uj. We 
find 

V($, A, a, uj, ...) = V($) + )-M 2 (<&)[\ 2 + a 2 + uj 2 ] + T($)\<ru + ... (5.42) 
As expected A, a, uj have the same effective mass 



M 2 ($) = _A e +2*/V6 

L 2 



and a unique trilinear coupling 



(5.43) 



(5.44) 



Along the CB flow 

M 2 (v) 



T(v) = 



-2/3 



(2-v) 



V2L 2 

The relevant wave equation for the inert scalars Q(x,v) = e ipx Q(v) reads 

.2 1 1 



(5.45) 



\v 1 — V J 



+ 



Q = 



(5.46) 



4£ 2 v 2 (l-v) v(l-v) 2 
with Q — A, a, u. The solution that vanishes at the singularity v — is 

Q(v) = v a (l - v)F(a, a; 2a; v) (5.47) 

where 10 

(5.48) 
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The inclusion of the fields A, a and u does not modify in any significant 
way the near boundary analysis performed in [1, 2]. In particular the renor- 
malized one-point function of the operator dual to any one of them (say A) 
is simply given by 

(O x (x)) = 2A (0) (x) (5.49) 

10 Using Pv°\z) = F(—u, v + 1; 1; (1 — z)/2), F(a, a; 2a; v) can be expressed in terms of 
associated Legendre functions, but little is gained in this rewriting. 
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Analytic continuation of F(a, a; 2a; v) reveals the non-local relation be- 
tween A( ) and A( ) at the linearized level and yields the two-point function 

N 2 

(O x (p)O x (q)) = 4— (2n) 4 5(p + q)iij(l) - ^{a)} (5.50) 

where the correct normalization of the action has been taken into account. 
The spectrum is continuous above a mass gap m = £/L 2 . Similar analysis 
leads to the bulk-to-boundary propagator (in momentum space) 

K p (v) = - v)F(a, a; 2a; v) (5.51) 

Taking the mixed derivative of (5.49) w.r.t. the sources of the remaining 
two fields (iT(o) and u>(o)) gives the desired three-point function 

(OM)OM)OM) = (2n)%Y.Pi) f Q dv^G(v)T(v)X{K P Xv) (5.52) 
Plugging in the relevant expressions yields 

r(aj) 2 r 1 

G CB (Pi^P2,P3) = «cB^EPi)II vto \ / dv ( 2 ~ v)\{v a 'F(a il a i - ) 2a i ]v) 

i i 1 \ la i) J i 

(5.53) 

where k cb = -(2N 2 /v / 3)(2tt) 2 (£ 4 /L 2 ). Since the SO (A) = SU(2) L x SU(2) R 
tensor structure is unique, Wigner-Eckart theorem gives 

(Q tLlR (pi)Q 3UR (p2)Q kLkR (P3)} = e t ^ kL e^ RkR g cB (p 1 ,p 2 ,p 3 ) (5.54) 

Although the integral cannot be performed elementarily, it is easy to see 
that G CB {pi,P2,Pz) is finite for generic momenta, as expected for operators 
of UV dimension A = 2, and vanishes at the threshold p 2 = —ml where 
dj = —1/2. Series integration yields a compact expression 

G CB (PliP2,P3) =«cbSII ' 



l + J2i(ai + ni) 2 + J2i(ai + rii) 

(5.55) 

where a S(J2iPi) is understood. 
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Direct comparison with the conformal three-point function of scalar pri- 
mary operators of dimension A = 2 computed in appendix A seems hopeless. 
However using the AdS / CFT correspondence for pure AdS we know that 

/dzd^x z 2 
~i5~n [ ;Z 2 + ( a ._ a ..)2]2 ( 5 - 56 ) 

can be written in momentum space, up to constants, as 

Qcft (Pi,P2,Ps) = K CFT 5(J2Pi) [ -j T (°)Il z2 ^o(zpi) (5.57) 

J z 

'I I 

where 

oo 2n 

M«) = E 7^7-K- 2 pog(«/2) - ^(n + 1)] (5.58) 

n=0 Z l 77 ^ 

is a modified Bessel function. Changing radial variable according to 11 

l-v = ^-\^ + 0^) (5.59) 

it is easy to convince oneself that the UV limit I — > 'explodes' the boundary 
region t; ~ 1 where indeed the asymptotics of K{p, v) coincides with the 
expansion of z 2 JC (zpi) and the trilinear coupling goes to its fixed point value 
T(0). This fixes the UV behavior including the overall normalization. 



6 Three-point function in the GPPZ flow 

Although the exact form of the fully non-polynomial potential of the M = 2 
gauged supergravity that governs the GPPZ flow is known 12 , it is not easy 
to isolate trilinear couplings that only involve inert scalars. Partly for this 
reason and partly to push our analysis another step forward, we are lead to 
consider the still relatively simple 3-point function 

G GPPZ (x 1 ,x 2 ,x 3 ) = {O a {x 1 )O a {x 2 )0 ( i ) {x z )) (6.60) 

11 The exact coordinate transformation with p = z 2 / 'L 2 can be found in [2]. 
12 We thank G. DalPAgata for disclosing its secret form to us. 
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where 0$ is the operator dual to the active scalar $ and O a is the operator 
dual to the inert scalar E. Q GPPZ vanishes at the superconformal fixed point 
in the UV by 5*0(6) R-symmetry arguments. 

The situation is subtler than in the previous Section. On the one hand 
$ mixes with the longitudinal and trace components of the metric. On the 
other hand a naive analysis inspired by the "dynamical scalar gauge" shows 
that the above 3-point function is logarithmically divergent as e — > if no new 
counterterm is generated when E is included. However, slightly at variance 
to what happens when only $ is present, we find two new logarithmically 
divergent counterterms. One is quartic in E the other is biquadratic in $ 
and E. 

In order to derive the precise form of these counterterms one has to per- 
form a detailed near boundary analysis of the bulk field equations. The 
results we need are 

$ = \^ - l R(0) ^) - ^(0) + 4))^ ( 6 - 61 ) 
/i(4) = -2(0(0)0(2) + o"(o)5"(2)) (6.62) 

9(4) = ^(9(2)9(2)) - 2(0(0)0(2) + <7( )<7(2)) + 4(0( O )0(2) + ^(0)^(2)) 

0(2) = -^(V( O )0(O) - (#(2))0(O)) - ^0(0) - 0(0)^0) (6.64) 

5"(2) = - J(V(o)<7 (0 ) - (g(2))(T(o)) + <rf ) - 0(o) ff (o) (6.65) 

(6.66) 

where R {0) = g^R^ij, gp) = g^gpyn, g(4) = g^g^a, h a) = 9(o) h m and 
{g(2)g(2)) = 9(o)9(2)jig[o)9(2)ki- 

Using Einstein equations in the convenient form 

iV = -2(dfld v $ + ^E^E + ^VG,J) (6.67) 

and including the Gibbons- Hawking boundary term, it is easy to see that the 
(regulated) on-shell action is simply given by 

S Icg = -\\ dpd 4 xVGV + ed e I ^Jg( x ,e) (6.68) 
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Inserting the data arising from the near boundary analysis gives 

f f 3 1 

SregbtO),^),^)^] = J d A X^/g^\-— + 7^(0(0) + 



^2 \ 



(6.69) 



+ loge 



1 

32 



R(0)ijR l (p) - 0-^(0) ) + 9 (0(O)V (o) 0(o) + 0-(O)V (o) <T( O )) 



+ O(e ) 



+ 7^(o)0(o) + ^(o)) - 2cr (0) + 2 ^(0) (7 (0) 
Once expressed in terms of the regulated fields, the counterterms read 

5 ct [ 7 ,$,E;c] = / p= / 4 ^{^-^W + ^($ 2 + S 2 ) (6.70) 
(W H $ + ^[ 7 ]$ 2 ) + 7 fEV [7] S + ii?[7]S 2 ) (6.71) 



-loge 



4 V 



4 V 



+2S 4 - 2£ 2 $ 2 + 1 (i^Wfrl - ^W 2 )] } (6.72) 

As anticipated the symmetry between <3> and £ is spoiled by the logarithmi- 
cally divergent quartic terms at the beginning of the last line. 
The one-point function of in the presence of sources is 



(O^) = [im{2(0 (2) + (2) ) - -0j o) - 41oge0 (O )<r 2 o) } 



(6.73) 



Despite appearance (O^) is finite as we will see. 

In order to compute the desired three-point function we need to determine 
the (non-local) dependence of 0( 2 ) on (X( ). Because of the mixing of the 
fluctuation ip of the active scalar $ with the longitudinal (if), radial (h rr ) 
and trace {h = g^hij) components of the metric one has to study the coupled 
field equations to linear order in <p, h rr , and H and to quadratic order in 
a. Following the steps detailed by [7, 8], but keeping h rr , one gets 



16 

h! H $V - 4A'h rr = J c 

3 

H" + AA'H' - -e- 2A {h + 2h„) = J H 
2 



(6.74) 
(6.75) 



3A'{ti - d 2 H) + -e- 2A d 2 h + AVhrr - 4$V - AV+tp = J G (6.76) 
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where d 2 = S^didj and 



Jc = ^(a'd t a) = (6.77) 



4e~ 2A / d i d j \ 4e~ 2A / d i d j 

J G = 2(a') 2 - 2e- 2A d i ad l a - 2M 2 a 2 = AT^ (6.79) 
Introducing the 'gauge invariant' combinations 

P = h+ M? R = K+2d \w) Q (6.80) 



one gets 

P' -4A'R= -J c (6.81) 
Q' + AA'Q - \e- 2A P - e~ 2A R = J H (6.82) 

3A'P' - 3A'd 2 Q + ^e~ 2A d 2 P + AVR = J G (6.83) 
Manipulating these equations as in [7, 1] eventually boils down to 

R'' + (2W H -4W)R' + (2W^W^-m 2 + ^-W 2 -^-WW H -e- 2A p 2 )R = J R 

(6.84) 

where the cr-dependent source 

J« = -4.V - - 8 , (--^ - 2.4' (-_) (6.85) 
turns out to be relatively simple thanks to the identity 

J' G +AA 'J G +^d 2 e- 2A J c -3A' 'J H +3A' ' J' c + 12(A') 2 J c = -2d r (M 2 )a 2 (6.86) 
that follows from 

V&Z^ = -\d„(M 2 )a 2 (6.87) 
where V# is the background covariant derivative. 
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Switching to the standard variable u and to momentum space yields 



R" + 



-R' + 



1 



2u -1 p 2 L T 



u(l — u) u(l — u) \u(l — u) 4 
The solution of the homogeneous equation is 



R = J R {u) (6. 



r(af)r(a^) to 



{n\y 



X 



x (1 - u) n [2^(ra + 1) - V(n + a?) - ^(n + aL w ) - log(l - «)] (6.89) 



where 



For u« 1, admits an expansion of the form 

R = i? (0 )(l - u) log(l -u) + R (0 )(l -«) + .. 

In the axial gauge /i rr = one finds 

4 4-1 

R (0) = -^0(2) + 0(2) -0(0)) #(0) = ~~yj^( 2 ) = ^| V ?O)0(O) + 

\/3 and 0( O ) 



To linear order around the GPPZ flow, where 0( O ) 
while ( 2) = cr (0) = 0, 

Differentiating i?( ) w.r.t. 0( O ), brings in an overall factor of p 2 that shows 
up in the two-point function 



(6.90) 
(6.91) 

(6.92) 
1/V3, 

(6.93) 



(0*(p)Ot(-p)) 



4^2 



p 



iP(af(p))+il>{a™( P ))-2i/>(l) 



(6.94) 



A similar though simpler analysis [7, 8] leads to the other two-point func- 



tion 



(O a (p)O a (-p)) = 



4^2 



P 



L 2 



^{at\p))+^\p)) 



(6.95) 
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where 

at\p) = \±\^-p 2 L 2 (6.96) 

Using 

00 z - 1 

^(z) = -7 + E 7 7T7 V ( 6 - 97 ) 

it is easy to identify the mass poles [11] and compute their positive (!) 
residues. 

For 0$, that belongs to the M = 1 'anomaly' multiplet A = E/Tr (Zf) 
dual to the active hypermultiplet, one has 

[m^Lf = A(k + l)(fc + 2) fc = 0,1,2,... . (6.98) 

For O a , that belongs to the M = 1 'Lagrangian' multiplet 5 = Tr {W 2 + ...) 
dual to the dilaton hypermultiplet, 

(m a L) 2 = Ak(k + 3) k = 0, 1, 2, . . . , (6.99) 

including a zero-mass pole. 

It is amusing to observe that the superglueball decay constants arising 
from the residues at the k th mass pole are the same for and a 

N 2 m 2 (k)(2k + 3) 
lM*OI = —^J} = IM*0I ! ( 6 - 10 °) 

We are ready to compute the three-point function. To quadratic order in 

<r's 



<^(P)) £ = ~R(o)(p) = du^G(u)K R (u,p)J R (u,p) (6.101) 

where, setting the active scalar to its background, 

J fl ( u ) = -16(</) 2 - 32(1 - u)a'a - 12a 2 (6.102) 

Differentiating (6.101) twice w.r.t. <7( ) one gets the logarithmically divergent 
integral 

rl-e 

G GPPZ (Pl,P2,P3) = K GPPZ 5C£2Pi)n T ( a +(Pi)) T ( a -(Pi)) J duU 3 (l-u) 2 X 

F R (p 3 ) [{F'MFM + F a ( Pl )F^p 2 )) - 2(1 - ujF'MF'M] (6.103) 
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rl-e 



where k gppz = (N 2 V^/2ti 2 L)(2tt) 4 and 

F R (p) = F(af(p),a W (p);3;u) (6.104) 

and 

F a (p) = F(a { °\p), a { °\p);3; u) (6.105) 

The logarithmic divergence of R(n\ is a contact term independent of the 
momenta. Its coefficient r] = 4v^3 is exactly cancelled by the last contact 
term in (0$) when 0( O ) is put to its background value of \/3. So that (0$) 
as well finite as they should. It is relatively easy to check 

that Q GPPZ vanish in the UV, i.e. at large momenta as expected on 5*0(6) 
R-symmetry grounds. 

Altough the above integral cannot be expressed in terms of elementary 
functions, it is not difficult to integrate the products of hypergeometric func- 
tions by series. Moreover the expressions drastically simplify when one goes 
on-shell i.e. p 2 = —m 2 and amputates the external legs thus effectively com- 
puting the irreducible vertex that enters the 3-body scattering amplitude or, 
equivalently, the superglueball decay amplitude. 

To this end observe that near the mass-shell where 

p 2 L 2 = -m(ki) 2 L 2 + a (6.106) 

a-(Pi) ~ — ki —r- r and a+(pi) ~ 3 + fa H — — — ^ (6.107) 

4(2^ + 3) +Kyj 4(2^ + 3) v ; 

This is true both for a and 0. Moreover, at the mass-shell the hypergeometric 

functions truncate and become degree k (Gegenbauer or Jacobi) polynomials. 

The mass poles of the connected 3-point function G GPPZ {jPi,P2,Pz) are 

exposed by the T's in front of the integral. Indeed 

r(a_(pO) - r (-ki - A J\ - ) - H* +1 4(2fr + 3) (6 _ 10g) 



4(2k t + 3) J kM 



while 



T(a + ( Pl )) « T (3 + h + 4(2A / + 3) j « (ki + 2)! (6.109) 

The relevant integrals are of the form 

f 1 du(l - M) 2 M Sn+1 [2n 1 n 2 (l - u) - (m + n 2 )u] = (6.110) 
Jo 

3\2n x n 2 2!(nx + n 2 ) 

(En + 2)(En + 3)(En + 4)(En + 5) ~ (£ n + 3)(E n + 4)(£ n + 5) 
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replacing the Pochhammer symbols in the hypergeometric functions with 
their explicit expressions 



{ ~ fa - th)\ {6 + k ^~ (k,i + 2)! (6 - m) 

and putting everything together yields 

A GPPZ (h, k 2 , k 3 ) = (2n ^ f[ \J (hi + ljfc + 2)(2fc, + 3) x (6.112) 

■_y«+m (^ + ni + 2 )!" 



NL 

i=i 



A / 6nin 2 -K+n 2 )(E» + 2) A 
n %\ (n s + 2)nU(En + a) l\ 



rii\(ni + l)\ (h -rii 



which has the expected mass dimension, [m] 1 and is correctly suppressed 
by the 1/JV factor at large N. Comparison with field theory results seems 
prohibitive for the time being since it would require a non-perturbative un- 
derstanding of the dynamics of J\f = 1 SYM which is well beyond reach. One 
may take the holographic result as a very precise prediction for the amplitude 
at strong 't Hooft coupling. 



7 Concluding remarks and speculations 

We have studied three-point functions of scalar operators along two holo- 
graphic RG flows. For the CB flow we have considered inert scalars and 
derived a compact expression that is finite and behaves correctly in the UV. 
In the GPPZ flow, though the initial situation was subtler because of the 
mixing of $ with the metric, we have been able to extract the irreducible 
vertices of three superglueballs. Thanks to the operator identity 

T> = ^ (7.113) 

with (5 = — y^, the very same modes couple to the trace of the stress tensor 
[1, 6]. This should elucidate the evolution with scale according to the holo- 
graphic renormalization group [34] . Extending the analysis to the transverse 
modes of the stress tensor or of the (broken) currents is straightforward in 
principle since the relevant trilinear bulk couplings are completely fixed by 
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gauge invariance. The expected form of these correlators should be 

(TMOMOM^ = M(£p*) J dw^G{w)Kl hkl {p u w) K7.114) 

i 

x[e~ 2A ^ w \p k 2 p 3 +P l 2 p k 3 ) + ^8 kl Ml{w)]K a (p 2 ,w)K a (p 3 ,w) + contact 

and 

(Ji(pi)O a (p 2 )0*(Ps))t = KjS(EPi) [ dwyjG^Kl^w) X 

i J 

(A -p k 3 )K a (p 2 ,w)K a (p 3 ,w) + contact (7.115) 

As expected, the integrals diverge but the contact terms suggested by holo- 
graphic renormalization cancel these divergences. Similar considerations ap- 
ply to three-point functions with more insertions of the stress tensor or of 
broken currents. As usual, bulk gauge invariance translates into Ward iden- 
tities of the boundary quantum field theory. The holographically computed 
one-point functions in the presence of sources do satisfy these Ward identities 
including anomalies and this carries over to higher point functions. 

It would be nice to address similar issues in fullfledged string solutions 
that display confinement and logarithmic corrections to the pure AdS behav- 
ior in the UV [26, 27, 28]. The results obtained by Krasnitz for the case of 
two-point functions in the KT flow are very encouraging in this respect [14]. 

As far as higher than three point correlation functions are concerned, if 
the miracle that in AdS transforms exchange diagrams into contact ones does 
not take place in non-trivial holographic RG flows then three-point functions 
will set the standard of a new holographic tradition, at least for a while. 
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Appendix: Conformal 3-point function in mo- 
mentum space 

In order to have a glimpse of the result of the integral for the holographic 
3-point correlation function in the CB flow and to eventually compare with 
not-so- standard field theory results it is convenient to compute the Fourier 
transform of the 3-point function of scalar primary operators of dimension 
A; i = 1,2,3. 

g cFT ( Xl ,x 2 ,x 3 ) = (0( Xl )0(x 2 )0(x 3 )) = ^ FT (x? 2 )- Al2 (^ 3 )- A23 (^3i)" A31 

(7.116) 

where x\ = (x t - xj) 2 and A^ = A; + Aj - E with E = (A 1 + A 2 + A 3 )/2. 
For simplicity we set the constant k cft = 1 henceforth. By definition 

S CFT (Pi,P2,p 3 ) = / ^x^d^Xse-^^+^^^ixuX^Xs) (7.117) 
In d dimensions, Feynman parametrization leads to the integral 
r , s s<sr ^ (s+d)/2 r(rf/2-A 23 )r(rf-E) 

G cft (puP2,ps) = ag» 2Sr (d/2)r(A 12 )r(A 13 ) x 

x f daF ^ - 1 a ^ p 1 - rll ~ a) T~T~\f^Z\ 

Jo 2 2 [(1 - a)(p 3 +p 2 ) 2 + (JP2\ 

where 

(1 - a)ap 2 3 

£.( a )=n \7 , \2~, 2 (7.118) 

(1 - a)(p 3 + p 2 ) 2 + api 

We have checked that this expression correctly factorizes in the extremal case 
Ai = A 2 + A 3 . 

Let us specialize to the case A; = 2 in d = 4, so that £ = 3 and A^- = 1. 
Putting 

ii 2 =p\lp\ li 3 = p\lp\ (7.119) 

one gets 

7T 7 / 2 

Q CFT (Pi,P2,Ps) = — —5{Y,Pi)A{V2,V3) (7.120) 
*Pi i 

where 

am ^ /- ^(i-'H-^+a-^i (7 . 121) 

jo (1 - a)a + a/i 2 + (1 - a)ii 3 
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is strikingly reminiscent of the box integral B(r,s) computed e.g. in [35]. 
Indeed A(/i 2 , A* 3 ) is tightly related to B(r, s). 

It is easy to check the symmetry A(/i 2 , A* 3 ) = A(fj, 3 , ji 2 ). It is slightly more 
laborious to express it in terms of logs and dilogs. To this end it is convenient 
to factorize the denominator as 

a(l -<t) + fi 2 a + (1 - 0> 3 = -(a - a + )(a - <r_) (7.122) 

where 



<7± = -[1 + /i 2 - /i 3 ± yj(l +fi 2 - Hs) 2 + 4/x 3 ] (7.123) 
Repeatedly and carefully (since they are valid for < x < 1) using 

Li 2 {x) + Li 2 (l -x) = - log(x) log(l + — (7.124) 



and 



where 13 



Lt 2 (x) + Li 2 = -\ log 2 (a:) (7.125) 



pi- dry 7" 

Li 2 {z) = ~ - log(l - za) = £ - 2 (7.126) 
Jo a n=1 n 2 

one finally gets the desired result 

a, \ 1 for- f 1 +^-^3-V^\ , or- + 
A{/jl 2: fi 3 ) = \2Li 2 I — ] + 2Li 2 



\ 2pL 2 J ~ \ 2pL 



- log(/J 2 ) log(/i 3 ) 
where 



log 



'fJ>2 + A*3 - 1 - V^" 



n 2 



(7.127) 



1/ = 1 + /is + A*3 - 2/^i - 2/i 2 - 2/i 2 /i 3 (7.128) 

A(fi 2 , /i 3 ) is manifestly symmetric and real for z/ > 0. Up to an overall 
constant, A(/i 2 , fj, 3 ) coincides with B(r,s) or better g(x±, x 2 , x 3 ) in [35] after 
replacing /j 2 and /j 3 with r and s, z.e. pi — > x 23 ,p 2 — > £ 31 , p 3 — > a?i 2 , 

13 In the notation of Abramowitz-Stegun [36] (pag.1004) Li^iz) = /(l — z). 
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Appendix B: Other useful formulae 

For the RG flows under consideration, homogeneous solution of the fluctua- 
tion equations may be expressed in terms of hypergeometric functions 

(f) p (w,x) = w a (l — w) l3 F(a, b; c; w) exp(ipx) (7.129) 

. We assume that z = represents the deep IR interior (singularity) and w = 
1— e the (regulated) UV boundary. In order to fix the overall normalization of 
the (regulated) bulk-to-boundary propagators K e (p,z) we need the analytic 
continuation of F(a, b, c; w) to w ~ 1. 
For c = a + b 

F(a , b;a + b; «,) = - g «(1 - «,)" x (7.130, 

x [log(l - w) - 2ip(n + 1) + ip(a + n) + ip(b + n)] 
For c = a + b — m 

, s r(m)T(a + 6 — m) (a — m) n (b — m) n 

x (1 - - (-rr(a + 6-m) ^ (a) n (b) n _ 

x [log(l — w) — ip(n + 1) — tfj(n + m + 1) + *0(a + n) + + n)\ 
For c = a + 6 + m 

j-t/ , MhM , T(m)T(a + b + m) ^ (a) n (b) n 

F(a, b;a + b + m ]W ) = na + ^ - ~ } E ~ "0 

(-)T(a + 6 + m) g (a + m) TO (6 + m) TO (i _ ^ )n+m ^ (? ^ 



r(a)r(6) ^ n!(n + m)! 
x [log(l — w) — ip(n + 1) — ip(n + m + 1) + ^(a + n + m) + ip(b + n + m)\ 

Correspondingly, the (regulated) bulk to boundary propagators read 



KM = _ ^ni-^rwr(6) 

loger(a + b) 
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K e (p,w) 




w a (l — w) /3 F(a + m)Y{b + m) 



F(a, b; a + b — m; w) 



(7.134) 



K £ (p,w) 



ePT(m)T(a + b + m) 



F(a, 6; a + b + m; ^J.135) 
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